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ABSTRACT 

This paper summarizes a three night observing campaign aimed at achieving 
milli-magnitude precision photometry of bright stars (V < 9.0) with the 1-meter 
Swope telescope at Las Campanas Observatory. The test targets were the main 
sequence stars HD205739 and HD135446. The results show that, by placing a 
concentric diaphragm in front of the aperture of the telescope, it is possible to 
avoid saturation and to achieve a photometric precision of 0.0008-0.0010 mag 
per data point with a cadence of less than 4 minutes. It is also possible to 
reach an overall precision of less that 0.0015 mags for time series of 6 hours or 
more. The photometric precision of this setup is only limited by scintillation. 
Scintillation could be reduced, and therefore the photometric precision could be 
further improved, by using a neutral density filter instead of the aperture stop. 
Given that the expected median depth of extrasolar planet transits of about 0.01 
mags, and their typical duration of several hours, the results of this paper show 
that 1-m telescopes equipped with standard CCDs can be used to detect planet 
transits as shallow as 0.002 mags around bright stars. 

Subject headings: techniques: photometric — stars: planetary systems. 

1. Introduction 

High precision photometry of bright stars is becoming a subject of increasing interest on 
extrasolar planet studies. Precision photometry complements the Doppler observations by 
establishing whether the observed radial velocity variations are caused by the reflex motion 
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induced by a planetary companion or by stellar magnetic activity (e.g., Henry et al. 2000a, 
Queloz et al. 2001, and Paulson et al., 2004). Precision photometry can also determine 
the presence or absence of transits in stars with known planets. The presence of a transit 
provides additional information about the planet by allowing the determination of its radius 
and its absolute mass (e.g., Charbonneau et al. 2000, Henry et al., 2000b, Sato et al. 2005, 
and Charbonneau et al. 2005), and also more precise orbital parameters of the system. 

Photometric precision of a few milli-magnitudes is routinely achieved using ground-based 
telescopes and standard CCDs for stars fainter than V = 9.0-10.0. Gilliland et al.(1991) 
showed that, using a 1-m telescope and a standard CCD, it is possible to achieve a precision 
of 0.0008 mags for 13th magnitude stars, in time series of about 9 hours with 2-minute 
candence. Everett & Howell (2001) have more recently demonstrated that it is possible 
to reach precisions of 0.00019 mags in a 4.5-hour time series for V = 14.0 stars with a 
0.9-m telescope and a mosaic CCD. However, for stars V = 9.0 mag or brighter milli-mag 
precision photometry basically remains the realm of photomultiplier tubes 1 . A typical CCD 
inmmediatly saturates when we place it in a 1-m class telescope and try to image a bright 
star. Saturation can be avoided by using a smaller telescope, but in that case atmospheric 
scintillation limits the photometric precision to 3-5 millimags. 

Efforts to reach milli-mag precision photometry of bright stars with CCDs on ground 
based 1-m class telescopes are now beginning to produce results comparable to those of 
photomultipliers tubes. For example, Bouchy et al. (2005) recently published the detection 
of an extrasolar planet transit around HD189733 with the 1.2 m telescope at the Haute- 
Provence Observatory in France, using a 1024 x 1024 SITe back-illuminated CCD. The depth 
of this transit is 0.03 mags and their photometric precision per data point varies between 
0.002 and 0.003 mags. Similarly, Charbonneau et al. (2006) have confirmed the 0.003 mag 
deep transit of HD149026b (Sato et al. 2005) using a 2-chip 2048 x 4608 CCD installed at 
the 48-inch (1.2m) telescope of the F. L. Whipple Observatory at Mount Hopkins, Arizona. 
Most efforts have been focused on the Northern Hemisphere, where small telescopes are more 
at hand. However, at least half of the sample of stars currently searched for planets lie at 
southern latitudes. In this paper I present the results of tests performed at the 40-inch (lm) 
telescope in Las Campanas Observatory, Chile. The tests were intended to achieve milli-mag 
precision photometry of southern bright currently being searched for planets by the radial 



1 The most successful ground-based experiment using photomultipliers is the automatic photometric tele- 
scopes (APTs) at Fairborn Observatory (Henry 1995a, Henry 1995b, Stammeir et al. 1997, Henry 1999, 
and Eaton, Henry & Fekel 2003). The APTs routinely achieve photometric precision of 0.001 mag for single 
observations of stars brighter than V = 8.5. They are also the only instruments that have maintained that 
precision for more than ten years, what makes them a superb tool. 
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velocity surveys. The selected targets were HD205739 (V = 8.5) and HD135446 (V = 8.2). 
Section 2 describes the instrumental setup and the observations. Section 3 summarizes the 
data analysis and presents the resulting light curves. Section 4 compares our results to the 
three currently known exoplanet transits around bright stars. Finally the summary of this 
work and future plans are given in section 5. 

2. Instrument Setup and Observations 

2.1. Instrument Setup 

The data was collected at the Henrietta Swope telescope, located in Las Campanas 
Observatory in Chile. The Swope is a 40-inch (lm) reflector with an f/7 Ritchey- Chretien 
optical design. It is currently equipped with a 2048 x 3150, 15 /iin pixel, SITe CCD that 
provides an unvignetted FoV of 14.8 x 22.8 arcminutes. The dynamic range of the CCD 
is 32727 ADU and it converts electrons into ADUs at a fixed gain of 2.5 e~/ADU. The 
readout time of the CCD on its 1 x 1 unbinned configuration is 128 sec. In addition, there 
is a minimum exposure time limit of 5 sec to ensure that the shutter has moved completely 
out of the way. 

With that setup stars brighter than V= 8.5 will unaviodably saturate in the 5 sec 
minimum exposure time. Saturation occurs even if we try to distribute the photons over a 
larger number of pixels by defocusing the telescope. 

One way to prevent saturation is to use a neutral density filter to block away a fraction 
of the photons arriving from the stars. Unfortunately no neutral density filter is currently 
available at Swope. Another way to reduce the number of photons impacting on the CCD 
is to place an aperture stop to reduce the photon collecting area of the telescope. 

With the help of Las Campanas Observatory personnel, I built a stop to reduce the 
effective aperture of the telescope from 0.589 m 2 to 0.146 m 2 . The stop, a 1 meter in 
diameter piece of plywood painted in black and with a central aperture of diameter 0.67 m, 
was placed in front of the aperture of the telescope as shown in figure 1. Objects appear 2.0 
magnitudes fainter with that setup, making it possible to image stars as bright as V = 6.5 
without saturating. 

To reduce in some amount the duty cycle of the observations, the readout of the CCD 
was limited to an area of 2048 x 2048 pixels (14.8 x 14.8 arcmin). That area is still large 
enough to simultaneously image the target star and at least one nearby comparison star of 
similar magnitude, but reduces the readout time to 90 sec, instead of the 128 sec that takes 
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to read out the full CCD. 



2.2. Observations 

The test targets for this work were the two southern stars HD205739 (V = 8.57; B-V = 
0.50; Sp. Type: F7V) and HD135446 (V = 8.20; B-V = 0.57, Sp. Type: G1.5V). Both stars 
have no previous record of any kind of intrinsic variability and are being searched for planets 
by the radial velocity surveys. HD205739 was observed during three nights, on UT 2005 Aug 
3-4, Aug 8-9, and Aug 9-10. HD135446 was observed on the nights of UT 2005 Aug 8-9 
and Aug 9-10. A log of the observations is presented in Table 1. HD135446 was monitored 
over the first three hours of each night. HD205739 was monitored over the second half of 
the nights (about 6 hours per night). The observations were all collected under photometric 
conditions, spanning over air masses between 1.01 and 1.63. 

The star used as comparison to HD205739 was HD205860, a star slightly brigther than 
HD205739 and of the same spectral type (V = 8.27; B-V = 0.50, Sp. Type: F7V). This last 
point is very fortunate, since it eliminates from the differential photometry the second order 
extinction effects introduced by differences in color between the comparison and the target 
stars. The separation between the two stars is Act = 46 s . 32 and A5 = - 10'. 31. In the case 
of HD135446 the comparison was HD135342 (V = 9.27; B-V = 0.51, Sp. Type: F6V). The 
separation between target and comparison is now Act = —25 s . 04, and AS = V.87. Both 
stars have also similar colors in this case, but the comparison is 1.1 magnitudes fainter than 
the target. 

The telescope was defocused so that the count level per pixel in both the target and 
comparison stars remained within the linearity limit of the CCD (< 25,000 ADU). The 
FWHM of the stars during the observations range between 15 and 19.5 pixels (6.5-8.5 arcsec) 
in the images of HD205739 and between 19 and 23.5 pixels (8.2-10.2 arcsec) in the images of 
HD 135446. While keeping the count levels below 25,000 ADU, the exposure times also had 
to be long enough to 1) ensure that enough photons were collected in both the targets and 
the comparison stars to keep the Poisson noise level below one part in a thousand, and 2) 
to reduce the errors introduced by scintillation. Typical exposure times were 20-22 seconds 
for HD205739 and 30-32 seconds for HD135446, respectively. 

Pixel to pixel sensitivity differences is one of the main precision limiting factors in CCD 
photometry. To minimize this contribution to the noise, the stars were placed over the same 
pixels each night, and then monitored with careful telescope guiding to ensure that the center 
of the stars was always located within the same ± 1-2 pixels. All the images were taken 
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using a V-band Johnson standard filter. 

Finally, several hundred bias and flat frames were collected each night to ensure an 
overall count level of over 10 6 photo-electrons per pixel in the final calibration frames. In 
this manner one can avoid introducing additional random noise in the calibrated images 
beyond the 1/1000 intended photometric precision. 

3. Analysis and Results 

All the images were bias substracted and flat-fielded using the same combined bias 
and combined flat frames. Those frames were generated by combining all the biases and 
flats collected in different nights, so that each of the final two calibration frames had an 
overall count level of at least 10 6 photons. That level of counts ensures that the random 
noise introduced during the calibration of the images, assumed to be Poisson noise, does not 
exceed 0.001 mags. 

The next step of the reduction processs consisted on performing aperture photometry 
of each target and their comparison star on the calibrated images. In order to optimize the 
results of the aperture photometry, I performed the photometry over series of apertures with 
radii varying between 12 and 35 pixels. The area used to compute the sky background around 
each star was always the same (an annulus of 15 pixels at a radius of 60 pixels centered on 
the stars). 

Differential photometry light curves were generated iteratively for all the possible com- 
binations of apertures of the target and the comparison star. The best light curve was then 
selected as the one resulting from the combination of apertures that gave the smallest value 
of the standard deviation (J ft/ fa where f t and f c are respectively the instrumental flux counts 
of the target and the comparison star. First, a prelimimary minimum of <7/t// c is sought by 
assigning aperture increments of 1 pixel per iteration. Once an approximate value of the 
minimum of <7/ t // c is found, one runs a second pass around that minimum to optimize the 
solution. In that second pass the aperture increments are 0.1 pixels instead of the 1.0 pixels 
in the first pass. The optimum apertures found were different for each target and each night. 
An example of the aperture versus Gft/fc plots generated to determine the best combination 
of apertures is illustrated in figure 2, where I show the results of the aperture photometry 
of HD205739 for the data collected on Aug 3-4. The optimum apertures in that case were 
14.8 pixels for the target and 15.6 pixels for the comparison star. The optimum apertures 
for each target in all the other nights are listed in Table 2. 

The final light curves are presented in figures 3 and 4. Figure 3 shows 16.95 hours of 
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photometric coverage of HD205739 over the three nights that this star was observed (Aug 
3-4, Aug 8-9, and Aug 9-10). Figure 4 shows the 6.2 hours of coverage of HD 135446 during 
Aug 8-9 and Aug 9-10. The data in both figures have been averaged into 2-point bins, 
resulting on a time of cadence of 3.5 minutes for HD205739 and 4.0 minutes in the case of 
HD135446. 

The average rms of the individual points range between 0.0008 and 0.0010 mags, de- 
pending on the night. The average standard deviation of the nightly light curves range 
between 0.00119 and 0.00130 mags for HD205739 and 0.00120 and 0.00162 for HD135446, 
indicating that both targets and their comparison stars were photometrically stable to at 
least those levels of precision during the nights that they were observed. 

The similarity in brightness and color of the targets to their comparison stars minimizes 
the appearance of differential extinction effects. These extinction effects are also reduced by 
limiting the photometric follow-up to airmasses smaller than 1.65. None of the light curves 
present differential extinction trends. I believe that the slight 0.003 mag slope between times 
2453586.80 and 2453586.87 HJD in the top light curve in figure 3 is real, after discarding 
extinction or background brightness gradient effects in the frames. Unfortunately, there are 
no follow-up observations to further investigate the cause of that slope. 

Scintillation is the limiting factor to the photometric precision of this setup. Equation 
10 from Dravis, Lindergren & Mezey (1998) yields a scintillation photometric variation of 
0.0012 mags for the average of two 25 sec exposures at an intermedian airmass of 1.3 2 . That 
value is in agreement with the standard deviation of the light curves in figures 3 and 4. 

4. Application to Extrasolar Planet Transits 

The expected median transit depth of giant close-in planets is 0.01 mags and the duration 
of those transits is typically of several hours. By comparing these numbers to the results 
presented in the previos section, one can conclude that a typical extrasolar planet transit 
would be easily detected around bright stars with the telescope setup and the observational 
plus data reduction strategies used Swope. 

To reinforce this statement, I provide in figure 5 an schematic representation of how 
the three known transits of extrasolar planets around brigth stars (V < 9) would appear 
when supperimposed on one of the light curves obtained in this work. The continuous line 



2 Assuming a telescope aperture diameter of 66 cm, an observatory height of 2100 m, and an atmospheric 
scale height of 8000 m. 
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simulates the transit of HD 149026b (Sato et al. 2005), the dashed line simulates the transit 
of HD209458b (Charbonneau et al. 2000), and the dotted line simulates the transit of 
HD 189733b (Bouchy et al. 2005). The transits have been represented by boxcar functions, 
instead of using a closer model to their real shape, to emphasize their duration At and their 
depth Amag as compared to the observed light curves, represented by the dots. The values 
of At and Amag have been estimated from the original papers and are summarized in Table 
3. 

The main conclusion from figure 5 is that the three extrasolar planet transits would 
have beeing easily detected by the setup that we are using at Swope. 

5. Summary and Conclusions 

This paper shows the results of a three-day test campaign aimed at achieving milli- 
magnitude precision photometry of bright southern stars (V < 9.0) with a 1-m telescope 
and a standard CCD. The targets for these tests were HD205739 and HD135446, two 8th 
magnitude stars with no previous record of photometric variability. 

Stars that bright would inevitably saturate given the collecting area of the telescope 
and the minimum exposure time of 5 sec imposed by the displacement rate of the CCD 
shutter. In the absence of neutral density filters at Swope, saturation can be avoided by 
implementing an aperture stop to reduce the aperture of the telescope, dimming the stars 
by two magnitudes, and making it possible to image objects as bright as V = 6.5 without 
saturating. 

The observations consisted of follow-up photometry of HD205739 over three nights, for 
a total of 16.95 hours. HD135446 was also monitored over two nights, for a total of 6.2 hours. 
Special care was taken to keep all sources of noise below the aimed photometric precision 
goal of 0.001 mags. The total number of counts in both the targets and the comparison stars 
was at least 10 6 photo-electrons. The level of counts per pixel never exceeded the linearity 
limit of the CCD. The level of counts in the combined bias and master flat-field calibration 
frames was also at least 10 6 photo-electrons. The stars were monitored only at airmasses 
under 1.65 to reduce differential extinction effects. The stars were also placed and kept 
over the same pixels each night to minimize errors caused by the pixel-to-pixel sensitivity 
variation of the CCD. Finally, we chose as comparison stars nearby objects with magnitudes 
and colors similar to the targets to minimize the effects of extinction. 

The aperture photometry was performed by iteratively generating differential photom- 
etry light curves for different combinations of apertures of the target and the comparison 
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star. The best light curves were then selected as the ones resulting from the combination of 
apertures that gave the smallest standard deviations. 

The resulting light curves show that it is possible to achieve precisions between 0.0008 
and 0.0010 mags in individual points with a cadence of less than 4 minutes. The results also 
show that it is possible to keep the level of precision of the light curves below 0.0013-0.0015 
for extended periods of time (at least 6 hours). 

The precision of the light curves is limited by atmospheric scintillation, which with the 
diaphramed aperture of the telescope and the current exposure times (20-30 sec) introduces 
an average photometric variation of 0.0012 mags. Scintillation depends heavily on the di- 
ameter of the telescope. Its effect can be therefore reduced to 0.0008 mags if we replace the 
aperture stop by an inexpensive neutral density filter. The neutral density filter will have 
the same advantage of dimming the stars by 2 magnitudes, as the aperture stop, but will 
not have the inconvenience of the lose in telescope aperture that enhances the atmospheric 
scintillation. 

These tests demonstrate that the Swope telescope can be used to obtain high precision 
photometry of southern bright stars, with very small alterations to its current setup. In 
particular Swope can be used to detect transits of close-in giant planets around bright planet- 
hosting stars. For example, Swope would had easily detected the transits of the planets 
HD189733b, with a depth of 0.030 mags, HD209458b, with a depth of 0.015 mags, and 
even the only 0.003-mag deep transit of HD149026b. There are currently no facilities in the 
Southern Hemisphere pursuing a targeted planet transit search among known planet-bearing 
stars. Swope is becoming therefore the first one of those facilities. 

Future plans include the replacement of the aperture stop by neutral density filters and 
the extension of these tests to other wavelengths, mainly UBVRI Johnson filter passbands. 
Another technical goal is to try to improve the duty cycle of the observations. In addition, 
we are starting a targeted search for transits around the southern stars with known planets 
that have not been monitored yet. 
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Table 1: Log of observations of the test target stars HD205739 and HD 135446 



Date(UT) Object Time(HJD-2453500) Airmass 



2005 Aug 3/4 HD205739 
2005 Aug 8/9 HD135446 

HD205739 
2005 Aug 9/10 HD135446 

HD205739 



86.62829-86.86426 
91.46557-91.60766 
91.63257-91.85669 
92.46533-92.58154 
92.62891-92.87488 



1.18-1.42 
1.01-1.63 

1.11- 1.51 
1.01-1.59 

1.12- 1.58 



Table 2: Optimum photometric apertures for each target in each night 



Target 


Date 


Optimum Apert 


Optimum Apert 


a f t /fJ 




(UT) 


of Target (pix) 


of Comparison (pix) 




HD 135446 


Aug 8/9 2005 


18.6 


18.8 


0.00506 




Aug 9/10 2005 


23.1 


22.1 


0.00501 


HD205738 


Aug 3/4 2005 


14.8 


15.6 


0.00405 




Aug 8/9 2005 


19.2 


18.9 


0.00452 




Aug 9/10 2005 


18.8 


18.4 


0.00235 



^The values of &ft/fc in this table include outlayers (points that deviate by more than 5a from the mean 
value of the light curve). Those outlayers were removed from the final light curves. 



Table 3: Duration (At) and depth (Amag) of the transits of HD149026b, HD209458b and 
HD189733b 



Planet 


At (hours) 


Amag (mags) 


HD149026b 


- 3.0 


0.003 


HD209458b 


- 3.0 


0.015 


HD189733b 


- 1.76 


0.030 



Fig. I. — Aperture stop to reduce the effective aperture of the Swope from 0.589 m 2 to 0.146 
m 2 . The stop is a 1 meter in diameter piece of plywood with a central aperture of diameter 
0.67 m. The stop reduces the brightness of the stars by two magnitudes, so it is possible to 
image stars as bright as V=6.5 without saturating. 
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Fig. 2. — Example of the aperture vs. Cft/fc plots to determine the best combination of 
photometric apertures for the target and comparison stars, as described in section 3. The 
plot shows the values of Gft/fc (y _ax is) f° r different apertures around the target HD205739 
(x-axis). Each line corresponds to a different aperture for the comparison star. The best 
photometry in this case is obtained when using an aperture of 14.8 pixels for the target 
and 15.6 pixels for the comparison star. The values of cr/t/fc m this plot include outlayers 
(points that deviate by more than 5<r from the mean of the light curve). Those outlayers 
were removed from the final light curves. 
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Fig. 3. — 16.95 hours of photometric coverage of HD205739 over three nights. The top, 
middle and bottom light curves corresopnd, respectively, to the data collected on UT 2005 
Aug 3-4, Aug 8-9, and Aug 9-10. The average rms of the individual points range between 
0.0008 and 0.0010 mags. The standard deviations of the nightly light curves are, from top to 
bottom, 0.00119, 0.00126, and 0.00130 magnitudes. The data have been binned into 2-points 
bins resulting on a time of cadence of 3.5 minutes. 
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Fig. 4. — 6.2 hour photometric coverage of HD135446 over two nights, Aug 8-9 (top) and 
Aug 9-10 (bottom). The average rms of the individual points range between 0.0008 and 
0.0010 mags. The standard deviations of the nightly light curves are respectively 0.00162 
(top) and 0.00120 magnitudes (bottom). The data have been binned into 2-points bins 
resulting on a time of cadence of 4.0 minutes. 
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Fig. 5. — Schematic representation of how the three known transits of extrasolar planets 
around bright stars (V < 9.0) would look like when superimposed on one of our light curves. 
The continuous line represents the transit of HD149026b (Sato et al. 2005) and the dashed 
and dotted lines represent respectviely the transits of HD209458b (Charbonneau et al. 2000) 
and HD189733b (Bouchy et al. 2005). The values of At and Amag for each transit are 
summarized in Table 3. 



